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. Abstract 

Constrained dynamics on finite-dimensional trial manifolds of quantum state 
vectors appears in time-dependent variational calculations. This work presents 
a requantization method, by which properties of the exact eigenstates can 
be retrieved from the phase portrait of such approximations. Applications 
^ ■ to the coherent states dynamics, small amplitude vibrations and collective 

CN ! rotations, show that this approach extends standard procedures of the quan- 



tum many-body theory, such as the method of the projection operators and 



the random phase approximation. 
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Introduction 

Constrained quantum dynamics usually takes into account self-consistent 
mean-fields in time- dependent variational calculations. Initially it was ex- 
pressed in the form of the time-dependent Hartree-Fock (TDHF) equations 
for many-fermion systems. However, the hybrid character of TDHF, includ- 
ing both quantal and classical aspects, leads to difficulties of interpretation 
pp. Thus, each TDHF solution is a quantum object, ensuring the optimal 
description of a time-dependent wave function within the manifold of Slater 
determinants. Moreover, the small amplitude TDHF vibrations around the 
critical points represented by solutions of the static Hartree-Fock (HF) equa- 
tions, correspond to the random phase approximation (RPA). Yet, important 
properties of the quantum dynamics are lost: the TDHF equations are non- 
linear, and tunneling effects are suppressed [2]. 

Such classical elements appear not only in TDHF. Any restriction to a 
manifold of trial states represents more than an approximation, because the 
result is partly outside the framework of quantum mechanics. Therefore, 
to get an effective approximation it is necessary to come back in the orig- 
inal Hilbert space by a "requantization" procedure. However, the classical 
aspects of the mean-field dynamics are useful for the separation of the col- 
lective variables [31 HI [5], [6] . 

The requantization problem is not yet completely solved, although im- 
portant partial results have been obtained. The methods of induced repre- 
sentations [7j or geometric quantization [S] can be used to treat globally the 
peculiar geometry of the collective models, but the constructed Hilbert space 
is physically independent of the initial one. A different procedure concerns 
the quantization of the closed orbits, by using a condition of periodicity and 
"gauge invariance" (GIPQ) [HI H0|- This formalism resemble the geomet- 
ric prequantization [TTJ, providing integrality constraints of Bohr-Wilson- 
Sommerfeld type. However, although GIPQ operates within the original 
Hilbert space, it is still far from complete requantization because the rela- 
tionship between the periodic trial states selected by such constraints, and 
the energy eigenstates, remains ambiguous. 

In this work approximate eigenstates associated to the quantized orbits 
of time-dependent trial functions are defined by the requantization proce- 
dure from [12]. Its relevance was tested on examples, by comparison with 
other methods and experimental data [13]. The first part contains a short 
description of the geometric structures related to the constrained quantum 
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dynamics. In the second part is presented the formula proposed for the 
requantization of the periodic orbits. On examples, it is shown that this 
formula embodies naturally both the RPA and the method of the projection 
operators. 

I. The constrained quantum dynamics 

The full quantum dynamics in the Hilbert space TL is given by the time- 
dependent Schrodinger equation (TDSE) 

ihd t \i>) = H\ip) (l) 

where H is the Hamiltonian operator, defined on a dense domain in TL. This 
equation can be related to the formalism of classical mechanics by its deriva- 
tion from the "variational principle" 8S\^\ = 0, with S^] = J dt{ip\iTid t — 
H\ip). However, it can also be seen as the flow 

\^)=X H ^) (2) 

of the quantum states G TL, under the action of the Hamiltonian field 
XhW = ( — i/h)H\4>) defined on TL by the relation ix H u H = dh. Here h 
denotes the Hamilton function h(ip) = (ip\H\ifj) , h : TL — > R, andu; w G fl 2 (Ti.) 
is the symplectic form ([H], 4.2): 

u H (X,Y) = 2hIm{X(iP)\Y(iP)) , X(ip),Y(ip) G T^H . (3) 

In the following 7i is presumed finite dimensional, although part of the results 
remain valid also in the infinite dimensional case. 

The invariance of h to the action on 7i of the group U(l), defined by 
& c \ip) = c\ip), c G U(l), leads to the conservation of the norm || \ip) \\ = 
yj Thus, it is possible to apply the Marsden-Weinstein theory to reduce 
the dynamics from TL to the projective space associated to TC ([15], 5.5.C). 

In terms of L = {\if)) G TL, = 1}, and of the projection n : L — > P H , 

n(|V» = [|V>>] = {e^), |V) e L, G [0, 2vr]} , (4) 

any \Z]) G L having as projection a trajectory "f t = [\Z^)] C of the 
reduced dynamical system, determines uniquely a solution of the TDSE in 
TC. This solution has the form 

\1> t ) = Q ef \Z?) = e?\Z?) , (5) 
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where / is given by 

f\Z?) = X S (Z?) - \%) . (6) 

Assuming further that H is independent of time, then h = (ip\H\ip) is a 
constant £ , and \ipt) from (J5]) becomes: 

|^) = e-s £t -/o ^1*1^1^) . (7) 

This result has a peculiar geometrical interpretation in the prequantization 
theory. Thus, by the smooth, free and proper action $ : U(l) x L — *■ L of 
the "gauge group" 2/ (1), L is a principal circle bundle over P n ([15], 4.1M). 
Each normalized state in Ji generates a one-dimensional Hilbert space, and 
by the natural action of the structure group on this space, L is also the unit 
section in the associated complex vector bundle [TB]. The 1-form a G 1 (L), 
a^(X) = —ih(ip\X), X G T^pL, (da = u) h \l), is a connection form on L, 
having as curvature the reduced symplectic form u G Q 2 (P-n), U*u = uj h \l- 
The covariant derivative [1 lj of a state \Z) G L along 7 C Pn, It = 

n(|^)),is 

^ = V,\Z)= l -a{\Z))\Z) = (Z\Z)\Z) . (8) 
Therefore, in ([7]) the state vector 

\~ Z ^ ) = e -!l^z7\d T \z7) ][Zl) (Q) 

is autoparallel along 7, while \ipt) corresponds to the lift of the Hamiltonian 
current from to the bundle L. 

The eigenstates of H are critical points of the reduced dynamics on P-h- 
In particular, the ground state \ip g ) can be found, up to a phase factor, as 
the minimum of h on L. This property is also used whenever \tp g ) is approx- 
imated by the constrained minimum |M) of h on a finite-dimensional "trial" 
manifold P C L. In static variational calculations such as HF or Hartree- 
Fock-Bogoliubov (HFB), the trial manifolds are orbits of the semisimple Lie 
groups SU(n) or SO(2n), respectively, generated by unitary action on high- 
est weight vectors for their representation in 7i [TTJ . Worth noting, the 
projection M = II(P) for an orbit of this type is a "coherent-state" Kahler 
submanifold of Pn [18J. 

If | M ) G P is an approximate ground-state, then also the constrained 
dynamics on P could be relevant for the description of the quantum sys- 
tem. Let P be a manifold of normalized trial functions P = {\Z%)} so that 
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M = n(P) is symplectic. The local coordinates of the point n(|Z^)) G M 
will be specified by 2N real parameters ). The Hamilton 

function on M, hu '■ M — > R is the restriction of h{ip) = (i/}\H\tp) to M. 
The symplectic form u M of M, U*uj m = ui n \p, defined by the restriction of 
uj n to P, can be expressed in the form 

2N 

uj m = ojijdx 1 A dx j , tOij = MmidiZldjZ) , (10) 
(<% = d/dx l ). Therefore, the Hamiltonian flow on M is given by the equations 



dhu 
dx 



The bundle structure II : L — > is reproduced at the level of the con- 
strained system only if M is quantizable, so that Lm = U(l)-P is a principal 
bundle over M. In this case the restriction of the 1-form ia/h G fi 1 (L) to L^/ 
yields in the local system of the trial functions |Z) G P the 1-form associated 
to the connection 9z G fi 1 (M), 

2N 

6 z = Y J {Z\d k \Z)dx k . (12) 

k=l 

In the prequantization theory, the connection provides the lift the Hamilto- 
nian flow from M to Lm- If %t is a solution of ( flTJ) . its lift pt|Z) is: 

p,|Z) = e-^ +ie '|^) , (13) 

where t 

8 = (Z$ t \H\Z St ) , ©t = / dr{Z iT \id T \Z iT ) . 

Jo 

The state pt|Z) represents a more appropriate approximation to an exact 
solution of TDSE than \Z t ) = \Z% t ). The same result was obtained previously 
0, [TU] by using the variational equation 

5 J C(x,x,r,r,<j),^)dt = (14) 

with C = (ty\ihdt — H\^) and = re l ^\Z,j. t ) . Therefore one can say that 
pt\Z) yields the best approximation of a TDSE solution in Lm- 
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It is interesting to remark that <d t coincides with the Berry's phase, de- 
rived first for the eigenfunctions of the adiabatic time- dependent Hamiltoni- 
ans [19l[20], then extended [21] and related to the quantization rules [22] . 

In general, the eigenstates of H are not elements of P, but nevertheless 
they influence the Hamiltonian flow on M. The correspondence between this 
flow and the energy eigenstates makes the object of requantization. 

The closed orbits on M can be quantized by using the condition of "gauge 
invariance"^ and periodicity (GIPQ) [9], [10] - This condition was obtained pre- 
suming that a periodic trial function resemble the most a time-independent 
eigenstate. However, beside the argument based on such similarity, GIPQ 
has a geometrical significance, revealed by the prequantization formalism 
[T5] . Thus, if C denotes the set of closed orbits j t = U(\Z)]) on M, then 
GIPQ approximates the stationary states by the autoparallel vectors: 

\Z t ) = e^\Zl) , 6* = f dr{Z1\id T \Z1) , 1 eC (15) 

Jo 

which are periodic. In the "gauge" of the functions \Z" f ) having the same 
period T 7 as 7, denoted |?7 7 ), the periodicity condition \U^) = \U^ +T ), 7 G C, 
leads to the "quantization rule" Gy 7 = 2irn, n G Z. Stated in terms of the 
1-form 8jj, this rule takes the form 

ijffcez. (16) 

Unlike the initial periodicity condition, this expression is gauge- dependent. 
Therefore, it is convenient to use the local relation oj m = —iTidOu to express 
( TT6l) in the intrinsic form 

7- f lu m G Z , 7 = da . (17) 

h J a 

Suppose now that the orbits of the Hamiltonian field on the energy surface 
Ps = hjj(£) are closed, and let uj^f be the restriction of u M to Pg. If ui^ 
is regular, then these orbits are also the leafs of the foliation T defined by 
the characteristic distribution of uo^ 1 [15]. In this case [23] the condition (fl6l) 
selects the energies £ for which the quotient spaces Ms = PejT are quanti- 
zable manifolds. 

2 The invariance up to a constant phase of \Zj) in ((9|) when is replaced by e lS ^ \Z]), 
with S(t) an arbitrary function of time. 
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II. The requantization of the periodic orbits 

If \Z t ) is a solution of some time-dependent mean- field equations, the eigen- 
values E n from H\ip n ) = E n \ip n ) can be approximated by £ n = (Z in \H\Z in ) 
for the orbits 7" G C selected by (1161) . However, if |Z 7 ") is taken as an 
approximation for \i/j n ), then its time-dependence, though periodic, leads to 
ambiguities in the evaluation of the transition amplitudes. To avoid such 
difficulties, in the following will be discussed the possibility of approximating 
\ip n ) by the time-averagqj |7 n ) of the periodic vector |Z 7 ™), 

m = ^idt\z^) (is) 

1 ™n J 



Example 1. The exact stationary solutions 

Let C = {7 = IT(|Z 7 )), \Z2) = e' iAt l % \Zl), t E R} be a regular orbit cylinder 
in P H [15] corresponding to the exact TDSE solutions. In this case each |Z 7 ) 
is quasiperiodic, and the related periodic gauge function is a "Bloch wave" 
in time, \U 1 ) = e %tlt ^ n \Z]). The quasienergy e 7 is defined up to integer 
multiples of Tiuj 1 = h/T 7 , and is chosen within the first "Brillouin zone" 
[— hu) J /2, hu 1 /2\. The condition selecting the quantized orbits 7 n becomes 
£ 7 = {Z2\H\Z^) = nhu~< + e 7 , n e Z, while 

| 7 ") = J- cf dte- it{ "- £ -< n)/h \Zf) (19) 

1 yTl J 

represents the projection of \Zq ) on the subspace spaned by the eigenstates 
with energy £ 7 n. The result is close to the ergodic mean theorem ([15J, 3.7.24) 
and can be easily illustrated for the harmonic oscillator, choosing the orbit 
cylinder associated to the Glauber coherent states. 

3 Knowing from [TU] that in some cases the time-averaged overlap between \Z 1 ) and 
\ip n ) attains a maximum for 7 = 7™, I used the time-average to estimate the order of 
magnitude of some matrix elements in a computer program [12] . However, it turned out 
that even the first significant digits were unexpectedly close to the experimental values, 
suggesting that the time average (|T8|) has its own relevance. The result was confirmed by 
an analytic model (M. Grigorescu, Rev. Roum. Phys. 34 1147 (1989)). In 1993 I found 
the reference S. Drozdz, M. Ploszajczak and E. Caurier, Ann. Phys. 171 108 (1986), in 
which a similar time-average was proposed. 
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Example 2. The harmonic oscillator 

For the one-dimensional harmonic oscillator H = h~uj(¥b + 1/2) with b = 
Jmuj /2%{x + ip/mu). Let P be the manifold of the Glauber coherent states 

P = {\Z) = e z&w6 |0> = e- |z|2/ V St |0), z G C, S|0) = 0} . (20) 

The trajectories on M = II(P) are given by z t = e~ luJt z , (Z t \id t \Z t ) = u>\z \ 2 , 
and the condition ( 1161) selects the values Zq , £ n : 

|4")| 2 =nGN , £ n = n u (\zi? ) \ 2 + ^)=hu(n+~) . (21) 

Because Zj = e~ ia3t ZQ , the period T of 7" is independent of n. Thus, the 
time-average of \Z? ) = exp(mu;t + z^b^ — z^*b)\0) expressed by 

|7»> = -f dt\Zf) = e -^j dte*»*ef* |0) (22) 

yields, up to a normalization factor, the eigenstate \ip n ) = (n!)~ 1//2 (6^) n |0). 
Worth noting, the "coherence" of the Glauber states makes the lift pt\Z) an 
exact solution of TDSE. 

Example 3. The Random Phase Approximation 

Let G be a semisimple Lie group with algebra g = T e G, and M = n(P) 
the Kahler submanifold of coherent states projected in from the group 
orbit P = {\Z) = U a \M),a G G, M G g*} of the highest weight vector 
|M) G H. In a local frame of 2n real parameters {£j,j = 1,2 AT}, (10) yields 
2uf k = ih{M\[VjU,V k U\M), and flTT} becomes 

2N 

ihY,£A M \lVft,V k U\M) = (M\{H u ,V k U\M) . (23) 
i=i 

Here "Dfc^ = U~ l dkU and P^/ = U~ l HU . In general, to find all periodic 
solutions of (1231 is a difficult task, but locally, around the minima of Hm — 
(M|P[/|M), the periodic orbits exist and can be obtained by integrating its 
linearized form. 

If A G g* denotes the set of roots for g, S = {m G A, (m, M) < 0}, and 
the operators {P m , Hm, m G A}, 

pt, = P_ m , PT m I M> = (m,M) I M> , 
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represent in Tt the Cartan-Weyl basis for g, then around a critical point 
Uo\M) G P the parameters £j can be chosen as the real and imaginary parts 
Xm, Vm of the complex variables z m , m G S defined by 

|Z) = £ Z |M) , £/, = C/ oe E meE (^^-^^- m ) . (24) 

With these parameters, the components of the form uf- in £/o|M) are = 
u yLy, = °> and 

< yn = -ft(M|p^tf,:D a _ n tf|M> = -M mn (M|[£ m ,£L n ]|M) . (25) 
Retaining only the terms linear in z, the equations of motion (I2"3"j) become 
^ m (M|[£ m ,£L m |M) = (M^iLjM) - (26) 
^(M\[[(z n E n -z:E^),H Uo ],E^]\M) . 

nGE 

The state \g) = Uq\M) is a critical point if the constant term in (|26|) vanishes, 

(M|[^ ,i_ m ]|M) = , VmGS , (27) 

which is the static equation used to find Uq. However, further will be taken 
Uq = I, presuming that if the minimum exists, then it represents both the 
ground state and |M) G P. Moreover, \g) is supposed to be invariant to the 
symmetry group Gs of H. Otherwise, the action of the identity component of 
Gs on \g) generates a whole critical manifold for Iim, related to the collective 
modes. Important examples are provided by the anisotropic or superfluid 
ground states. These are not invariant to the action of the rotation group 
.50(3, R), respectively of the "gauge" group U{1) generated by the particle 
number operator, and the collective behavior shows up in rotational bands. 

Let |M) be a Gs-invariant minimum, and consider the closed orbits as- 
sociated to the normal vibration modes. These are specified by the complex 
amplitudes {X m ,F m ,m G S} and the oscillation period T, of the general 
periodic solution 

zl = X m e~^ + Y m e wt , m G S, u — — . (28) 

Replacing this expression in (|26l) . a time-independent system is obtained: 

(M| [[#,£+] -TwoB\E m ]\U) =0 , Vm G A, (m, M) ^ , (29) 
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where ffl denotes the sum 

# = £ (X m E m - Y*E_ m ) . (30) 

mGS 

For (PED, 6»[/ = (M|£H<9 r £>|M) from (H5J is dominated by the term 

^2(M\[E m ,E_ m }\M)Im(z m z* m ) = (31) 

mGS 

-^^(M|[E m ,E_ m ]|M)(|X m | 2 -|y m | 2 ) , 

mGS 

and the quantization condition for the amplitudes X m , Y m becomes 

- ^(M|[ J E m)J E_ m ]|M)(|X m | 2 -|F m | 2 )=n, n = l,2,... (32) 

mGS 

The result makes sense only if the quantized amplitudes are within the range 
of the linear approximation. Assuming this to be true for n = 1, the au- 
toparallel section with period T becomes 

\Z?)=e iut e e - iutm - eiut *\M) (33) 
and the approximate stationary state given by the time-average is 

|w) = i^>* «#|M) • (34) 

In the HF case G = SU(n), and the set of operators associated to the roots 
mGS contains A(n — A) operators {c\cj,A + 1 < i < n, 1 < j < A} 
expressed in terms of the single-particle creation and annihilation fermion 
operators {c|,q,z = l,n}. In this representation fl29l) becomes the RPA 
equation in the particle-hole channel, and for the first excited state (n = 1) 
the quantization condition fl32l) reduces to the normalization equation for the 
RPA operator^ (M\[B,&]\M) = 1. 

Similarly, by choosing G = SO(2n) the RPA equations in the particle- 
particle channeO are obtained. 

4 This relation between quantization and normalization was obtained using a different 
method by S. Levit, J. W. Negele and Z. Paltiel, Phys. Rev. C 21 1603 (1980), Eq. 3.39. 

5 If |M) is the particle vacuum, for m 6 S there are n(n — 1) operators {cfcj, i, j — l,n}. 
An example is presented in M. Grigorescu and E. Iancu, Z. Phys. A 337 139 (1990). 



10 



Example 4- The requantization of the periodic orbits in action- angle coordi- 
nates 

Let the dynamical system denned on M = n(P) be completely integrable 
in the action-angle coordinates (1,0), y5 = (if 1 , <p 2 , ip N ) G T^, I = 
(Ji, I2, In) G R ■ The manifold P will be represented as 

P = {\Z(0I)) = e-i^=^\Z(O,I))} (35) 

with [J^ J k ] = 0. If there is no degeneracy the closed orbits correspond to 
the fundamental cycles {7&,fc = l,iV} of M, 0*\-y k = Sj k X k . The condition 
(flE]) leads to the quantized expectation values 

(Z^\j k \Z lk ) =mh ,m G Z (36) 

and further to a quantized set of action variables I. 
The autoparallel vectors defined by (|T5|) are 

I^Tkj = e ^ fc (t) e -^ fc (*)4| Z 7 fc (o,/)) (37) 

and their average over the orbit 7™: 

| 7 -) = _L Idipe-^ Jk - mh) \Z^{0,I)} (38) 

is the projection of the state |Z 7 ™(0, /)) on the subspace spaned by the 
eigenstates of J k with the eigenvalue mh. 

In applications, constraining manifolds P of the form fl35l) appear when 
the variational equation 5(Z\H\Z) = on a general trial manifold V has a 
"deformed" (symmetry-breaking) solution \g). Let be the maximal torus 
of the symmetry group G$ of H. The ground-state \g) G V is presumed 
non- invariant to the action of any subgroup of T N , and is used to define the 
point \Z(0,0)) = |#) of P C V. The states \Z(<p,I)) are related to |Z(0,/)> 
by the action of T , and therefore to define P is necessary to specify the 
leafs parameterized by the action variables. Because 

2uy 4 = 2hIm(d^Z(0I)\d Ik Z(0I)) = d Ik (Z(0,T)\Jj\Z(0,I)) , 

the canonical expression 2uJ v jj k = 5j k is obtained if 

I k = (Z(<p,I)\J k \Z(<p,I)) = (Z(0, 1)1^(0, 1)) . 
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Moreover, because [H, J k \ = the Hamiltonian 



h M = (Z{(p,I)\H\Z{<pJ)) = (Z(0,I)\H\Z(0,I)) 

depends only on /. Thus, for a given J, the state \Z(0, 1)) can be taken as 
the minimum of hu = (Z\H\Z) on the sub manifold of V selected by the 
set of N conditions Ik = (Z\Jk\Z), k — 1,N. This constrained variational 
problem reduces to a normal one for the modified Hamiltonian 

N 

H- X = H-J2 A fc 4 (39) 
k=i 

where A = {A fc , k — 1, N} are the Lagrange multipliers. Denoting by \g x ) the 
solution of the variational equation 5(Z\H~ X \Z) = in V, and by A fc (J) the 
functions of / determined by the system of implicit equations Ik = (9\\Jk\9\), 
k — l,N, the result is \Z(0, 1)) = \g~ x(}) ). 

Discussion and conclusions 

The natural restriction of the quantum dynamics to the manifold of normal- 
ized states, related to the phase invariance, was used as a model also for the 
finite-dimensional, quantizable, trial manifolds. 

The correspondence between the constrained dynamics and the exact so- 
lutions of TDSE was studied for Hamiltonian operators independent of time. 
It was shown that the closed orbits on the trial manifold can be related to 
the stationary states by a procedure extending the prequantization formalism 
and the GIPQ method. This procedure requires: 

• The selection of the "quantizable" closed orbits lying under periodic 
autoparallel sections. 

• The calculation of approximate stationary states by the time-average 
of the selected autoparallel sections. 

In applications the choice of the constraining manifolds plays a central 
role. Particularly important appear to be the exact orbit cylinders and the 
manifolds associated to the symmetry-breaking ground states. Depending 
on this choice, the proposed procedure appears as a fine instrument either 
for the study of the exact quantum dynamics, or of the collective modes. 
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In particular, it shows the common nature of different techniques such as 
the restoration of broken symmetries by using projection operators, and the 
random phase approximation. 

It is important to remark that for the harmonic oscillator, among all trial 
manifolds the Glauber coherent states are distinguished because the con- 
strained dynamics coincides with the one of the classical system. Therefore, 
in this case the problem of quantization can be stated as the requantization 
of the constrained dynamics, rather than its correspondence with an isomor- 
phic, but different quantum framework. 

These results indicate the consistency and physical relevance of the pro- 
posed procedure. However, this should be considered only as a part of a 
general requantization formalism. A straightforward extension concerns the 
requantization of the orbits dense on invariant phase-space tori, and will be 
considered in a subsequent paper. 
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